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Abstract 

We discuss application of the physical QCD effective charge ay, defined 
via the heavy-quark potential, in perturbative calculations at next-to-leading 
order. When coupled with the Brodsky-Lepage-Mackenzie prescription for 
fixing the renormalization scales, the resulting series are automatically and 
naturally scale and scheme independent, and represent unambiguous predic- 
tions of perturbative QCD. We consider in detail such commensurate scale 
relations for the e + e~ annihilation ratio R e + e - and the Bjorken sum rule. In 
both cases the improved predictions are in excellent agreement with experi- 
ment. 
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I. INTRODUCTION 



One of the most important problems in making reliable predictions in perturbative QCD 
is dealing with the dependence of the truncated perturbative series on the choice of renor- 
malization scale \x and scheme s for the QCD coupling a s (/i). Consider a physical quantity 
O, computed in perturbation theory and truncated at next-to-leading order (NLO) in a s : 



l + (AM+BMn f )^- + - 

7T 



(1-1) 



where n/ is the effective number of quark flavors. The finite-order expression depends on 
both \i and the choice of scheme used to define the coupling. In fact, Eqn. fll.lj) can be made 
to take on essentially any value by varying \x and the renormalization scheme, which are a 
•priori completely arbitrary. The scale/scheme problem is that of choosing /i and the scheme 
s in an "optimal" way, so that an unambiguous theoretical prediction, ideally including some 
plausible estimate of theoretical uncertainties, can be made.[] 

For any given observable there is no rigorously correct way to make this choice in gen- 
eral. However, a particular prescription may be supported to a greater or lesser degree 
by general theoretical arguments and, a posteriori, by its success in practical applications. 
From these perspectives, a particularly successful method for choosing the renormalization 
scale is that proposed by Brodsky, Lepage and MacKenzie [3J. In the BLM procedure, the 
renormalization scales are chosen such that all vacuum polarization effects from fermion 
loops are absorbed into the running couplings. A principal motivation for this choice is that 
it reduces to the correct prescription in the case of Abelian gauge theory. Furthermore, the 
BLM scales are physical in the sense that they typically reflect the mean virtuality of the 
gluon propagators. Another important advantage of the method is that it "pre-sums" the 
large and strongly divergent terms in the pQCD series which grow as n\(a s (3o) n , i.e., the 
infrared renormalons associated with coupling constant renormalization. 

Dependence on the renormalization scheme can be avoided by considering relations be- 
tween physical observables only. By the general principles of renormalization theory, such 
a relation must be independent of any theoretical conventions, in particular the choice of 
scheme in the definition of a s . A relation between physical quantities in which the BLM 
method has been used to fix the renormalization scales is known as a "commensurate scale 
relation" (CSR) [0]. An important example is the generalized Crewther relation 0|§, in 
which the radiative corrections to the Bjorken sum rule for deep inelastic lepton-proton 
scattering at a given momentum transfer Q are predicted from measurements of the e + e~ 
annihilation cross section at a commensurate energy scale y/s oc Q. 



The precise meaning of "optimal" in this context is connected to the minimization of remainders 
for the truncated series. As is well known, perturbation series in QCD are asymptotic, and thus 
there is an optimum number of terms that should be computed for a given observable. In general, 
very little is known about the remainders in pQCD; however, if we assume that pQCD series are 
sign-alternating, then the remainder can be estimated by the first neglected (or last included) term. 
This term can take on essentially any value, however, by simply varying the scale and scheme, and 
thus its minimization is meaningless without invoking additional criteria. 
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A useful tool in these analyses is the concept of an "effective charge." Any perturba- 
tively calculable observable can be used to define an effective charge by incorporating the 
entire radiative correction into its definition. Since such a charge is itself a physical observ- 
able, perturbation theory in terms of it, with the BLM prescription setting the scales, is 
automatically renormalization scale- and scheme- independent. 

In this paper we shall use the heavy quark potential to define an effective QCD coupling 
ay-, and construct scale-commensurate expansions of various other QCD observables in 
terms of it. A recent calculation of the heavy quark potential at NNLO j| allows the 
relevant BLM scales to be determined through NLO. The resulting relations can be tested 
directly for agreement with available data, and in addition may be used to study various 
phenomenological forms for the heavy quark potential at moderate to low Q 2 . 

We begin by outlining the BLM approach and the idea of commensurate scale relations. 
We also introduce physical effective charges asociated with the heavy quark potential, the 
e + e~ annihilation cross section and the Bjorken sum rule. In section III we then construct 
the NLO scale-commensurate expansions of these observables in terms of ay, and compare 
the results to the available data using a simple parameterization for ay which is fit to a 
lattice calculation. In general the agreement is excellent. In section IV we present some 
discussion of the results and our conclusions. 



II. BLM SCALE FIXING 



At lowest order the BLM approach is straightforward to motivate. The term involving 
rif in Eqn. ( |1.1|) arises solely from quark loops in vacuum polarization diagrams. In QED 
these are the only contributions responsible for the running of the coupling, and thus it is 
natural to absorb them into the definition of the coupling. The BLM procedure is the analog 
of this approach in QCD. Specifically, we rewrite Eqn. ( |1 . 1|) in the form 



O = a s (ii) 



1 - 



'3/3 -Bi(/i)\ a s (/i) 



+ 



7T 
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(2.1) 



correct to order a 2 , where /3q = 11 — 2n//3 is the lowest-order QCD beta function. The first 
term in square brackets can then be absorbed by a redefinition of the renormalization scale 
in the leading-order coupling, using 



a s (fi*) = a s (fi) 



1 ^— — In (// Ifi) + 



(2.2) 



That is, the BLM procedure consists of defining the prediction for O at this order to be 

335i(»\ a s (n* 



O = a 8 (ji* 



A x (ii) 



+ 



7T 



where 



li = lie 



3-Bi(,u) 



(2.3) 



(2.4) 



Note that knowledge of the NLO term in the expansion is necessary to fix the scale at LO. 
Thus the scale occurring in the highest term in the expansion will in general be unknown. 
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A natural prescription is to set this scale to be the same as that in the next-to-highest-order 
term. 

A very important feature of this prescription is that fi* is actually independent of \i. 
(This follows from considering the \i dependence of -£>i(/i). For a detailed discussion of this 
point, see Ref. ][]].) Thus pQCD predictions using the BLM procedure are unambiguous. 

The same basic idea can be extended to higher orders, by systematically shifting rif 
dependence into the renormalization scales order by order. Full details of this procedure 
may be found in Refs. |5|U. The result is that a general expansion 

+ {Al + Bl n,) + (A 2 + B 2 n, + C 2 „J) fcMf +■■■ (2.5) 

is replaced by a series of the form 

^ + Al {^y + A 2 {^flJ + .... (2.6) 

In general a different scale appears at each order in perturbation theory, and the BLM 
scales themselves are power series in the coupling a s . In addition, the coefficients A n are 
independent of n/ (by construction), and so the form of the expansion is unchanged as 
momenta vary across quark mass thresholds. All effects due to quark loops in vacuum 
polarization diagrams are automatically incorporated into the effective couplings. 

As discussed above, one motivation for this prescription is that it reduces to the cor- 
rect result in the case of QED. In addition, when combined with the idea of commensurate 
scale relations, the BLM method can be shown to be consistent with the generalized renor- 
malization group invariance of Stiickelberg and Peterman J7J, in which one considers "flow 
equations" both in /x and in the parameters that define the scheme ||. This is not necessarily 
true of other methods for determining the scales. 

A very natural way of implementing the CSR idea is to introduce a physical effective 
charge, defined via some convenient observable, for use as an expansion parameter. An 
expansion of a physical quantity in terms of such a charge is a relation between observables 
and therefore must be independent of theoretical conventions, such as the renormalization 
scheme, to any fixed order of perturbation theory. A particularly useful scheme is furnished 
by the heavy quark potential V(Q 2 ), which can be identified as the two-particle-irreducible 
amplitude for the scattering of an infinitely heavy quark and antiquark at momentum trans- 
fer t = —Q 2 . The relation 

V{Q2) = _ AirC F ay(Q) m 

with Cf = (N 2 — l)/2N c = 4/3, then defines the effective charge ay(Q)- This coupling 
provides a physically-based alternative to the usual MS scheme. The other physical charges 
we shall consider here are ocr, defined via the total e + e _ — > X cross section: 

i?(^3Ee;(l + ^) , (2.8) 

and a gi , defined by the radiative correction to the Bjorken Sum Rule: 
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dx 



g r(x,Q 2 )-gr(x,Q< 



9a 



gv 



1 - 



a ai(Q) 
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(2.9) 



The perturbative expansions for these quantities through NNLO may be found in Refs. || 
and PJlOll, respectively. 

Such physical couplings are of course renormalization-group-invariant, i.e., fid ay jdfi — 0. 
However, the dependence of ay(Q) on Q is controlled by an equation which is formally 
identical to the usual RG equation. Since ay is dimensionless we must have 



ay = ay 



'Q 



Then ftday jdfi = implies 







Q-^a v (Q) = (3 s (a & 



a s {fi) 



day 
da„ 



Pv(<*v) 



where 
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(2.10) 



(2.11) 



(2.12) 



This is formally a change of scheme, so that the first two coefficients f3y =11 — 2n//3 and 



(3y^ = 102 — 38rif/3 in the perturbative expansion of (3y are the standard ones. 



III. QCD PERTURBATION THEORY AND ay 

A. BLM Scale Fixing for ay 

The calculation of the heavy quark potential at NNLO in Ref. 0] allows the BLM pro- 
cedure to be applied through NLO in commensurate scale relations involving ay. As a first 
step, we may apply the BLM procedure to fix the renormalization scales in the expression 
for ay in terms of the conventional MS coupling. The result is 

<*v{Q) _ ajrsiQy) + Ay f <*m(Qv) \ 2 , By f *Ms(Qv*) Y , . . . ) (31) 



where 



7T 7T \ 7T / \ 7T 



A v = -hj A (3.2) 

Bv = m _ ii C , + _ LA d + (-™ + L\A CaCf (3.3) 

v V144 4 S3 8 64 / A V 192 4 S V v ; 

ln(g* v /Q) = -| (3.4) 
o 

217 21 /105 9 \C F 

ln(Q v /Q) = -— + — Cs + tttt - -Ca 7T" ( 3 - 5 ) 



192 16 V128 8 / C A 
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and Ca = N c . As discussed above, we take Qy * = Qy at this order. 

It is also useful to invert this, and express aj^ itself in terms of ay. In this case we 
obtain 

^ = 2^ + ^(2^) 2 +Bw (M|ii)y + ..., ( 3. 6) 

where 

Ams = \ca (3.7) 



HQjzg/Q) = \ (3-9) 
ln(Q^/Q) = 192 + 16^ 3 + ( 128 ~~ 8^) C~ A ' (3 ' 10) 



B. e + e Annihilation Cross Section 

We next present the NNLO scale-commensurate expansion of cxr in terms of ay. This 
is obtained by applying the BLM procedure at NLO to the expansion of each of these 
observables in the MS scheme, and then algebraically eliminating ajjg. The result is 



<*r(Q) _ MQr) , Ar ( MQr) Y , Br 



7T 



7T 



7T 



+ 



7T 



(3-11) 



where (for N c = 3) 
Ar 



25 
12 



97 
72 



9 



10 



11 



-7T 



Ca 



64^ + d(r) \1U 6 



v e 2 



HQr/Q) 



23 

-- + 2C 3 + (33-2n / ) 

157 233 „ 

+ — Cs - 2C 5 • 



119 7T 2 7 2 2 

864 + 72 ~ <T 3 + 3^ 3 



60 50 



'tty(Q) 



7T 



(3-12) 

(3.13) 

(3-14) 
(3.15) 



In Eqn. (3.13), d(r) is the dimension of the quark representation, i.e., 3 for SU(3). This 
relation represents an unambiguous, fundamental test of perturbative QCD which is inde- 
pendent of renormalization scale or scheme. 
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FIG. 1. The momentum dependence of n>f(Q z 



In order to make a comparison to experimental data, we will introduce a parameterization 
of ay which is fit to lattice data in the moderate- to high-Q 2 regime. Specifically, we 
take 



ocv{Q) 



An 



(3.16) 



Asymptotically this reproduces the perturbative coupling, while the effective "gluon mass" 
m g results in ay becoming essentially frozen for Q 2 < 4m 2 ,. This form can be motivated on 
various theoretical grounds [jTJ] , and it has also been successful in phenomenological analyses 
0- 



The parameters Ay and m 2 have been determined in Ref. 
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by fitting to a lattice 
using 



calculation of V(Q 2 ) [Tl|| at relatively high Q 2 and to a value of a^ advocated in [|n 
Eqn. ( fmp at LO. They were found to be Ay = 0.16 GeV and m 2 = 0.19 GeV 2 . 

Note that in the beta function (5q we use a "smeared" function for the number of flavors, 
although this only affects the low-energy regime where several quark flavor thresholds occur. 
This function is 



n f (Q 2 



E 



dz- 



QQ 2 z 2 (l-z) 
i) + Q 2 z{\- 



(3.17) 



and is motivated in Ref. 15 



The integration over z in Eqn. (3.17) leads to the explicit 
representation] 2 ] of the function that is identical to the Q 2 logarithmic derivative of the one- 
loop massive /^-function presented in Ref. |15|]. In Fig. p] we show rif(Q 2 ) in the low-energy 



2 Note that Ji dz 6 ^^ 



1-Q X + , ln / (l+^/Hl 

1 UX ^ (l+4x)V2 m \(l+4x) 1 / 2 -l 
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FIG. 2. The effective charge ay, as given by Eqn. (3.16) 

region. We have taken m s = 0.15 GeV, m c = 1.9 GeV, mb = 4.5 GeV for the quark masses. 
The resulting ay is shown in Fig. |[ 

Note also that for low Q 2 the couplings, although frozen, are large. Thus the NLO 
and higher-order terms in the CSRs are large, and they do not give accurate results at 
low scales. In addition, higher-twist contributions to the effective charges, which are not 
reflected in CSRs relating them, may be expected to be important for low Q 2 . However, 
series expansions in terms of physical charges are likely to be more convergent than those 
cast in terms of unphysical couplings such as cxj^g, which is singular at finite scales.^ Thus 
it is quite possible that expansions of the type we are considering can be extended to lower 
physical scales than series written in terms of ajjg. In any case, we will not be directly 
concerned with the low-Q 2 regime here. 

Before discussing the results, it is useful to understand what improvements we can expect 
from the commensurate scale relations. First of all, of course, we have a scale-independent 
result, so aesthetically we have an advantage over the conventional treatment. Moreover, 
because of this we expect our result to be numerically more accurate than previous results 
with the scale fixed to certain value. The main applicability and usefulness of commensurate 
scale relations is for the intermediate energy regime. Pertubation theory is valid only above 
the characteristic QCD scale Aqcd, and since the commensurate scale analysis crucially 
depends on the validity of perturbation theory, we don't expect much improvement in the 
very low energy regime. Furthermore, in the high energy limit the residual scale dependent 
terms go to zero, so scale relations are meaningless. The e + e~ annihilation data, as well 
as the Bjorken sum rule data presented in the next section, lies in the intermediate energy 
regime where we expect improved predictions. 



3 For example, in the 't Hooft scheme ctjfq has a simple pole at Q = A-tfo. 
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Two additional modifications of Eqn. ( |3.11|) were performed before comparing with 
data. First, we have included the leading-order electroweak corrections to account for the 
Z° current, which is particularly important above 30 GeV. In addition we have included the 
charm and bottom mass corrections, which are important in the range 3-15 GeV. The effect 
of these modifications is to replace the factor J2 q e 2 , in Eqn. (|2.8|) by 



\ 



Am 2 



1 + 



2m 2 



+ 2Re(r)e,44 1 + 



2m 2 



+ \r\ 2 ((4) 2 + {c\f) (4) 



(3.18) 




4m 2 ' 



Q 2 



where 



r 



C A 



I 3 — 2e q sm 2 9 w 



2sin 2 6'w — - 



1 

2 ' 



y^GM 2 



Q*-M 2 +iM z T z \e 2 ) 

q_ 

sm 2 26 w (Q 2 - M 2 Z + iM z T z ) ' 



(3.19) 



Here I3 is the third component of the weak isospin of the quark coupled to Z° and the weak 
mixing angle 9w is given by sin 2 9w = 0.22. The mass and the decay width of Z° are given 
by M z = 91.2 GeV and T z = 2.5 GeV, respectively. 

In Fig. |3|, we show the commensurate scale result (|3.11|) along with a representative 
subset of the available data [[II]] in the energy range 8-60 GeV. We find our results to be 
in excellent agreement with the data, as well as the standard QCD predictions quoted by 
the Particle Data Group [17] with the scale fixed to a certain value (A^jg = 0.25GeV). In 
Fig. H, we show our theoretical prediction and the data in the 2-7.5 GeV range. Again, 
we find very nice agreement with the data, particularly considering that we have neglected 
corrections from the J/ip(lS), ip(2S), and other vector meson resonances. Note that the 
data for 3.6 GeV < Q < 7.5 GeV has been subtracted by 0.84^ 



l l-4m 2 /Q 2 (l + 2m 2 T /Q 2 ) 
to account for hadronic production that proceeds via tau lepton pairs, which the early 
experiments did not distinguish from quark-hadron processes. The factor 0.84 = 1 — (2/5) 2 
is the probability that either tau will decay to hadrons. 



C. Bjorken Sum Rule 



Finally we present the scale-commensurate expansion of the Bjorken sum rule in terms 
of ay at NNLO. The result is 
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FIG. 3. The scale-commensurate expansion of an in terms of ay in the high energy regime. 
The solid line is given by Eqn. (3.11); the dashed line is the prediction quoted by the PDG; the 
dotted line is the leading order result (with mass and electroweak corrections). 
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FIG. 5. The dotted line shows the leading-order prediction for the Bjorken sum rule while the 
solid line includes the scale-commensurate corrections through NNLO. Also shown are data from 
Ref. [18] (crosses), from the E154 Collaboration [19] (circles), and from the SMC Collaboration 
[20] (stars). 



a gi {Q) _ a v (Q* g 
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where 




+ Ag [-^S-L\ +B g ^ZJLM + ... ; ( 3. 20) 



An = — (3.21) 
9 12 v ' 

131 27 7 9 , . 

B Q = 7i 2 + —7i* 3.22 

9 72 8 64 K J 

1Q1 Fi 

HQT/Q) = - + ^C 5 • (3.24) 

In Fig. ^| we show the commensurate scale result to NNLO and the leading order pertur- 
bative result with the five currently available data points. This plot strongly suggests that 
the higher order pQCD corrections do indeed give the correct convergence to the physical 
result. Our results may also be compared with an analysis of the Bjorken sum rule |T^ 



using so-called analytic perturbation theory (APT) [21j . In Ref. [0 , the authors show that 



by requiring the QCD coupling a s to be analytic, thereby removing unphysical singularities, 
they can obtain approximately scheme independent results. Their plot of the correction to 
the Bjorken sum rule, a gi /ir, is very similar to what we obtain using commensurate scale 
relations. 
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IV. CONCLUSIONS 



In this paper, we have applied the physical QCD effective charge ay, defined by the heavy 
quark potential, in calculations of the e + e~ annihilation cross section and the Bjorken sum 
rule. Following the BLM procedure, we derived the NNLO scale-commensurate expansions 
of cxr and a 9l in terms of ay and used these expansions to numerically compute the e + e~ 
annihilation cross section and the Bjorken sum rule. Using a phenomenological form for 
the effective charge ay [Eqn. ( |3.16| )1 which is consistent with the lattice determination 
of the heavy quark potential, we obtain excellent agreement between our results and the 
experimental data in both cases. Furthermore, because of the scale independence, we trust 
that our results are numerically more accurate than previous results with the scale fixed to 
a certain value. The application of scale-commensurate expansions to other observables is 
forthcoming. 
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